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Magnetic flux detection with an Andreev Quantum Dot 
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The charge of the subgap states in an Andreev quantum dot (AQD; this is a quantum dot inserted into a 
superconducting loop) is very sensitive to the magnetic flux threading the loop. We study the sensitivity of 
this device as a function of its parameters for the limit of a large superconducting gap A. In our analysis, we 
account for the effects of a weak Coulomb interaction within the dot. We discuss the suitability of this setup 
as a device detecting weak magnetic fields. 
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Introduction. The Joscphson effect [1] has been in- 
tensively studied during the past 45 years; its main char- 
acteristic is the presence of a tunable non-dissipative 
current when two bulk superconductors are joined via a 
normal or insulating layer and subjected to a supercon- 
ducting phase difference tp. Recently, it has been real- 
ized that in a metallic junction the charge of the normal 
island in between the superconducting leads depends on 
the superconducting phase difference tp as well [2,3]. 
This dependence is sufficiently strong [3] to use this ef- 
fect in a magnetic flux detector, although our estimates 
below give a sensitivity somewhat below the sensitivity 
of the best SQUIDs. 

Usually, small magnetic fields are measured 
by superconducting quantum interference devices 
(SQUIDs) [4, 5]. While SQUIDs are based on the 
dependence of the Josephson current on the supercon- 
ducting phase difference tp (and hence on the magnetic 
flux $ threading the loop), here we propose to use 
|the charge-dependence in an Andreev quantum dot 
for the flux measurement. As shown in Ref. [3], the 
charge Q of a single-channel Andreev quantum dot can 
be fractional — |e| < Q < \e\ and depends on tp (here 
e = — |e| is the charge of one electron). 

The charge of an Andreev quantum dot can be 
measured by a sensitive charge detector, e.g., by a 
single-electron transistor (SET). Today, the best sin- 
gle electron transistors have a sensitivity of the or- 
der of 10~ 5 |e|/VHz (e.g., see [6]). Using results of 
Ref. [3] , simple estimates tell that an AQD can convert 
a change in flux (5$ to a change in charge SQ with a ra- 
tio 6Q/8& — 2|e|/$0j where $ = 2-jvh/2\e\ is the super- 
conducting flux. Assuming a superconducting loop area 
~ 1 mm 2 , we obtain the sensitivity 1(U 14 T/vTIz, which 
is comparable with the sensitivity 10~ 14 H- 10 -15 T/ vTIz 



of today's best SQUIDs [4,5]. Below, we study in detail 
the sensitivity ratio 5Q/S&. 

Setup. Our Andreev quantum dot is realized by 
a small metallic dot connecting two superconducting 
banks joined in a loop, see Fig. [TJ Our AQD is as- 
sumed to be a quasi one-dimensional normal metal (N) 
island separated from the superconductors (S) by thin 
insulator layers (I) , generating normal scattering on top 
of the Andreev scattering characteristic of the normal- 
metal superconductor junction. The position of the nor- 
mal resonance in this SINIS system can be tuned by 
the gate voltage V g applied to the normal region of the 
AQD. The magnetic flux threading the loop $ induces 
a superconducting phase drop tp across the AQD. Since 

Superconducting 
loop 




1 ' sadovsky@itp.ac.ru 



to SET 



Fig. \T\ Andreev quantum dot inserted into the super- 
conducting loop. The Andreev quantum dot is con- 
nected to a single electron transistor (SET) and a gate 
electrode through capacitive couplings. The flux $ pro- 
duces a phase difference tp = 27r"l>/$o across the An- 
dreev quantum dot. The charge of the AQD can be 
tuned by the gate voltage V g and the flux $ threading 
the loop. 
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the phase drop in the bulk superconductor is negligible 
as compared to the phase drop ip across the AQD one 
may relate the latter to the flux $ threading the loop, 
<p = In order to measure the charge trapped 

on the AQD, a single electron transistor is capacitively 
coupled to the normal metal island. Experimentally, 
such AQDs have recently been fabricated by coupling 
carbon nanotubes to superconducting banks [7-10]. In 
the following, we concentrate on the properties of the 
key element in the setup — the Andreev quantum dot. 

Energy and charge of the AQD without Coulomb 
interaction The Andreev states give rise to new op- 
portunities for tunable Josephson devices, e.g., the 
Josephson transistor [11-13]; here, we are interested in 
their charge properties. We will consider the case of one 
transverse channel such that the problem effectively be- 
comes one dimensional. We consider the case of a large 
separation S N between the resonances in the associated 
NININ problem (where the superconductors S have been 
replaced by normal metal leads N), <5 N 3> A, such that a 
single Andreev level £ A is trapped within the gap region. 
We are interested in sufficiently well isolated dots with 
a small width T N of the associated NININ resonance, 
T N <C A. In this section, we neglect charging effects 
E c = 0. In summary, our device operates with energy 
scales r„ « A « i B . 

The resonances in the NININ setup derive from the 
eigenvalue problem 7i * = E^> with H = -ft 2 <9 2 /2™ + 
U{x) — e F with the potential U{x) — U ps ^(x + L/2) + 
U ps , 2 {x — L/2)] + eV g 9(L/2 — \x\)] describing two point- 
scatterers 2 ) (with transmission and reflection ampli- 
tudes T ; V V', i? ; 1/2 e xr ; R l = l - T u I = 1, 2) and 
the effect of the gate potential V g , which we assume to 
be small as compared to the particle's energy E (mea- 
sured from the band bottom in the leads), eV g <C E. 
Resonances then appear at energies E n — El(jvk — 
xl/2 - X2/2) 2 ; they are separated by <5™ = (E n+1 - 
E n -i)/2 2E n /n and are characterized by the width 
r N = TS^/ny/k, where e L = h 2 /2mL 2 . The bias V g 
shifts the resonances by eV g ; we denote the position of 
the n-th resonance relative to e F by e™ = E n + eV g — e F . 
In the following, we choose a specific resonance in the 
gap by selecting an appropriate n and drop the index n, 

t N ' t N , (J N > N , 1 N f L N . 

We go from a normal- to an Andreev dot by replacing 
the normal leads with superconducting ones. In order 
to include Andreev scattering in the SINIS setup, we 
have to solve the Bogoliubov-de Gennes equations (we 
choose states with e A ^ 0) 

2 )The Heaviside function 6(x) = 1 for x > and 6(x) = for 
x < 0. 
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with the pairing potential A(x) = A[0(— x — 
L/2)e-^/ 2 + 9(x - L/2)e^/ 2 ]; u(x) and v(x) are 
the electron- and hole-like components of the wave 
function. The discrete states trapped below the gap 
derive from the quantization condition (in Andreev 
approximation) 

(Ri +R 2 )cos(2ny^ -4y / R 1 R 2 sin 2 a cos^ttj-^) + 
+T 1 T 2 cosip = cos^2a-27r^-^+i?ii? 2 cos(2a+27r^) . 

The phase a = arccos(e A /A) is acquired at an ideal 
NS boundary due to Andreev reflection with ip = 0; 
the above formula can be directly obtained using results 
from Refs. [13,14]. 

We concentrate on the regime T N , |e N | -c A, the so- 
called A — > 00 limit. In this limit, the quantization 
condition can be expanded and we obtain the expres- 
sion (A is the asymmetry parameter) 



£a — -\/e 2 + £ 2 , 



where 



£2 
2 



^ + A 2 , 
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A = 



2s/T^T 2 
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The energy e A of the Andreev state is phase sensitive 
when e N is close to the chemical potential, |e N | < T N , 
which can be achieved by tuning the gate potential V g . 
In the limit A — > 00, both the u(x) and v(x) components 
of the wave function are nonzero only in the normal re- 
gion, 



u(x) 
v(x) 



0. 



-ik c x 



\x\ > L/2, 
\x\ < L/2, 



where fc Ci h = [2m(e F ± e?/)] 1 / 2 /h are the wave vectors 
of electrons and holes, respectively. The coefficients are 
defined by C^ h = C^ h = [(1 ± e N /e A )/4L] 1 /2. 

The ground state of the system is the state |0) with 
energy 

£0 — £n - £a (5) 

(counted from the Fermi energy e F ), where we have sub- 
tracted the energy of filled resonances below the Fermi 
surface; the latter are not modified by the supercon- 
ductivity in the leads and hence do not depend on the 
phase (p. The first excited state with one Bogoliubov 
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quasiparticle is doubly degenerate in spin = aJ|0), 
= at |0) and has energy e\ = £q + £ a = %• The dou- 



form the entire basis of the system's Hilbert space 3 - 1 . 
The interaction is given by the operator 



bly excited state with two quasiparticles |2) = a|a||0) 
has an energy e 2 — £o + 2e A = £ N + £ A - 

The charge of the state \u) {v = 0, If, 2) 
can be obtained by differentiation of the correspond- 
ing energy £„ with respect to the gate voltage, q v = 
de u /dV g , or by averaging the charge operator Q — 

e J2<r I-l/2 ^t( x )^<r( x ) dx over the state \v), q v = 
(v\Q\v). Both methods give the identical results 



V = E c 



I. 

,2 ' 



(7) 



'('-*)• 



q_2 = e 



Given the basis with these four states, we can diago- 
nalize the Hamiltonian exactly. The non-zero matrix 
elements of the operator V are 

V m = E c {ql+ql 2 )/e\ V n = E c , 

V22 = E c (ql + q 2 02 )/e 2 , V 02 = 2E cqt) 2/e. (8) 

( 1 + — ] . (6) The energy levels are defined by the eigenvalue problem 



Below, we will also need the off-diagonal matrix ele- 
ments of the charge operator Q; the only non- vanishing 
term is q 2 = (0|Q|2) = e(l - e 2 N /el) 1 / 2 . 



AQD with Coulomb interaction In order to find 
the effect of weak Coulomb interaction E c <C A in the 
limit r N , |e N | <C A, we can disregard the continuous 
states with energies above the superconducting gap A 
and assume that the four levels of the discrete spectrum 
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Fig. [2] Energies E- (solid line), E\ (dashed line), and 
E+ (dotted line) are plotted versus the position of nor- 
mal resonance. All energies are given in units of e r , cf. 
|[3j). The Coulomb energy is E c = for (a), E c = £r 
for (b), E c — 2e r for (c), and E c = 3e r for (d). In ac- 
cordance with formula the doublet region appears 
when E c > £r, see (b-d). In the filled region the ground 
state of the system is a doublet; the width of this region 
is 2{E C — e^.) 1 / 2 , the edges of this region are spread due 
to the finite temperature 0. 
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where £„=£„ + V vu , 



0, It 



0, 



li, 2. The energy of 



the level with one Bogoliubov quasiparticle |1) is given 
by the (shifted) constant 



Ei — e„ 



E r 



(9) 



and does not mix with the other states; furthermore, the 
spin degeneracy of this Kramers doublet remains. The 
ground state |0) and the doubly excited state \2) mix 
due to Coulomb interaction and produce two new states, 



the singlet states |— ) and \-\ 
Dt/Df = -V 02 /(e ~ E±), 
energies of these new states are 



|±) = D*|0>+i?ri2), 
1. The 



\Df\ 2 



E±=e N + 2E C ± ^(e N + 2E c ) 2 +e 2 



(10) 



The energies of the doublet and singlet states de- 
pend on £ N and ip in a different way and may cross; 
thus the ground state can be formed by either the sin- 
glet |— ) or by the doublet The state |+) always 
remains the second excited state, see Fig. [21 When 
E c > £r ^ T N A/2 = E* (with A the asymmetry param- 
eter) the ground state is the doublet 1 1) in the region 



- 2E C - ^E 2 - 4 < £n < -2E C + y/E% - e 2 T (11) 

and remains |— ) at all other values of £ N [15]. 

The origin of this level crossing can be traced to the 
different shifts in energies with E c : While E\ is shifted 
up by E c , E— quickly approaches with increasing E c ■ 



3 'ln realistic nanodevices the Coulomb energy can be larger 
then r N and smaller or of the order of 5 N , but in principle can 
be made much smaller than both <5 N and A (see the discussion 
in [3,7]). 
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Superconducting phase, (p 



e N + 2E C =2.0E' C 
£ N +2E c = l.5E' c 

e N +2£ c = 1.0£c 
S N + 2E C = 0.5E' C 
E N +2E r =0 



e N + 2E c = -0.5Ec 
8 N +2£ , C = -1.0£ C * 
8 N +2£' C = -1.5£ C * 
e N +2£ = -2.0£ c * 



Fig. [3] Equilibrium charge Q cq (|15p versus supercon- 
ducting phase difference tp. In (a) and (b) the temper- 
ature is zero (i.e., Q eq represents ground state charge), 
in (c) and (d) the temperature is Q = Q.IEq, where 
E* c = V N A/2. The Coulomb energy is E c = l.2E* c for 
(a) and (c), E c = 2.QE* C for (b) and (d). The asym- 
metry level of the dot is A = 0.2. The features in the 
center of the plots corresponds to the Kramers doublet 
region pip. In (c) and (d) the border of the doublet 
region is smoothed by the temperature G. 



Note that the terms oc q 2 2 and oc (702 in the matrix ele- 
ments Vuv> lead to the crossing of the energies E- and 
Ei , while preventing the crossing of the level E + with 
the others. 

At the edge of the region (fTTj) a sharp singlet to dou- 
blet crossover takes place, with a jump appearing as a 
function of e N (V g ) or e r (<p) in the charge of the Andreev 
dot and in the current across (see below). The charges 
of the new states (// = 1, ±) can be calculated as in 
the previous section, = dE^/dVg, and are given by 



i± 



2E C , 



v/(e N + 2E c y + e\ 



Qi 



(12) 



The charge Q\ is integer and does not fluctuate; the 
charges Q± are fractional in the region |e N + 2E C \ ~ e r 
and fluctuate strongly (see also the discussion of fluc- 



tuations in Ref. [3] where Coulomb effects have been 
ignored) 



SQ ± ^[(±Q 2 \±)~(±Q\±) 2 ] 1 / 2 ^ 



v/(e N + 2E c y + e\ 



(13) 



Note that the Coulomb interaction merely shifts the 
regime of e N where the charges Q± are fractional. Ev- 
erywhere outside the doublet region the ground state 
charge is given by Q-, while within the Kramers dou- 
blet region the charge is pinned to the value Qi = e. As 
illustrated in Figs. [3K and [3b, for E c > E* a sharp 
crossover occurs and the charge jumps by the value 
SQcr = Q- — Qi- This jump is smeared at finite tem- 
peratures, see Figs. [3]; and[3ji. The groundstate charge 
is 

E N 



Q 



ss -e-e-^9[E c <E A ], 

E A 



where E A = [(e N + 2E c ) 2 +S 2 ] 1 / 2 and £ N = e N + 2E c de- 
notes the energy of the shifted normal state resonance. 
The equilibrium charge at finite temperature is 



Qeq — 



Q _ e -E-/e + 2 Q ie -*Vc- 



l+e 



-E+/Q 



(15) 



e -E-/e + 2e--Ei/© + e ~ B +/ e 

here and below we set Boltzmann's constant k B = 1. 
The equilibrium charge as a function of the supercon- 
ducting phase tp — 27r<!>/$o is shown in Fig. [3] 

The currents in the states are defined by rela- 
tionship J fj, = dE^/dQ which provides the results 



J± = T 



2ir T 2 sin tp 



$ 16E A 
The groundstate current is 



Ji = 0. 



(16) 



2ir T 2 sin tp 
$0" 16E A 



>[E C <E A }; 



note that the current vanishes throughout the doublet 
region. The thermal equilibrium current is 



Jot 



j_ e -E-/e + J+e - 



E+/e 



e -E-/e + 2e --Bi/e + e -E+/e ■ 



(18) 



Differential sensitivity The differential sensitivity 
of the equilibrium charge to the magnetic flux thread- 
ing the superconducting loop is defined by the absolute 
value of the derivative dQ oq / <9$ taken at the given value 
of flux, 4 ) S = \dQ eq /d$\. By using (I5J) we obtain 



S = 



F dQ dF e 



(19) 

4 'Note that the sensitivity of the charge-to-flux convertor S = 
Sfr^Q coincides with the voltage- to-current sensitivity of the 
Josephson transistor described in Ref. [13] Sv— >J = \9Jcq/9V g \. 
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where Q = (Q + — Q_)/2, the derivative 
2tt E N Tlsm(f 



the function 



dQ 



-E+/B 



16£? 



(20) 



-B_/e 



e -B_/e + 2e -Bi/e 



e --E+/0 " 

sinh(£J A /6) 
'cosh(^ A /e) + e £ c/e^ 



and its derivative 
<9$ 



sinh(£; A /e) + 1 
[cosh(£ A /e) + e £ c/e]2 



J_ 



(21) 



(22) 



As illustrated in Fig. [3] there are two intervals where 
the Qcqiv) dependence is steep. As ip increases from 
<p = 0, the charge increases (decreases) and reaches a 
maximum (minimum). For E c < E* the maximum 
(minimum) of the charge is always at if = 7T, while 
for E c > E* the extremum splits and a second inter- 
val with a steep dependence Q e q(cp) emerges in between 
the two new extrema. The first interval (interval I in 
what follows) corresponds to the singlet state of the 
AQD, the second (interval II in what follows) corre- 
sponds to the doublet state. We start with a descrip- 
tion of the first interval. We fix the parameters T N , 
A, and E c and search for the maximum sensitivity S 
as a function of e N and ip. The non-trivial symmetries 
Q eq (ip, e N ) = Q cq (2Tr-ip, e N ), Q cq (p, e N )-Q eq (cp, 0) = 
— Q eq (tp, — e N — 4:E C ) + Q cq (ip, 0) allow us to restrict the 
search to the region < ip ^ 7T, e N + 2E C ^ 0. Subse- 
quently, we analyze the maximum as a function of E c 
keeping A and T N constant. 

Interval I: For E c < [3(1 + A 2 ) /{I + 2A 2 )} 1 / 2 E* and 
zero temperature = the sensitivity is determined by 
the derivative dQ/d<$> O- The function \dQ/d<S>\ has a 
maximum at e N + 2E c = [(1 + A 2 ) /(I + 2A 2 )] 1 ^ 2 E* and 
tp = tt — 2 arcsinL4/ (1 + 2A 2 ) 1 / 2 ], where the differential 
sensitivity is given by 

^m^v e 



2tt 



$o 6V3AV1 + A 2 ' 



(23) 



One observes that the smaller A is, the larger is the sen- 
sitivity. In other words, a symmetric SINIS structure 
provides a better sensitivity S^ X (A — > 0) — > oo, but at 
the same time the region in <p with this large sensitivity 
vanishes as A — > 0. When <C E* the sensitivity is 
nearly independent of temperature. 

In the opposite case -Be > [3(1+A 2 )/(1+2A 2 )] 1 / 2 E* 
the doublet region covers all of the interval I and the 
maximum at zero temperature is always realized at the 



tN 

e° 

x 0.8 
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| 0.4 











© - 0.2£ c * 






- QAE* 






& = 0.65c* 






Interval II 












~ - _ Interval I 









^ 0.2 0.4 0.6 0.8 1 

Coulomb energy, i? c /r N 

Fig. [4] Maximum of the differential sensitivity S ma x 
(absolute value) in the interval I (dashed lines) and in 
the interval II (solid lines) versus Coulomb energy E c 
at the asymmetry level A = 0.2, E C /T N =0.1. The 
temperature varies from O = 0.2E C up to 6 = 0.6E C . 

edge of the doublet region (|ll|l . with a sensitivity given 
by 



S" 1 = 



,2tt 



$0 48\/3£ r 



x sj 2(A 2 - A + 1)3/ 2 - (A + 1)(A - 2)(2A - 1) (24) 

realized at e N + 2E C = (r N /2){[2A - 1 + (A 2 - A + 
l) 1/2 ]/3} 1/2 and p = 2arccos{[A + 1 - (A 2 - A + 

2 )/(r N /2) 2 . This re- 



l)l/2]/ 3 }l/2 ) where A = (#2 

suit reduces to 

ol 



E* 



,27r_r|_ 
1 $ 16-B 2 



(25) 



in the limit E c 3> T N , and remains approximately cor- 
rect for E c T N /2. For E c 3> r N , the maximum sen- 
sitivity is reached at e N + 2E C « E c — r 2 /16i? c and 

^«7r/2 + r 2 /i6£; 2 . 

Interval II: At zero temperature there is a jump in 
the charge at the edges of interval II and thus the sen- 
sitivity diverges in these points. A finite temperature 
smears the jump and the sensitivity becomes finite. If 
E c 3> O, T N , E*, the sensitivity S reaches the maxi- 



mum near the point £r, 
equals to 



all 

w max 



2E C — E c . ip — 7r/2 where it 



,2tt r 



(26) 

The expression for S}^ ax is too cumbersome for an ar- 
bitrary Coulomb energy E c and we plot the numerical 
result S^ ax (E c ) in Fig. [4] In the same plot, we also 
present the maxima of the sensitivity from the interval 
I. One easily notes that for a large Coulomb interaction 
the charge jump smeared by temperature provides the 
sharper Q oq {p) dependence. 
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Conclusion. In this article, we have pointed out that 
the ^-dependence of the charge trapped within an An- 
dreev quantum dot may be used for the implementation 
of a new type of magnetometer which operates along the 
pathway 'magnetic flux-AQD charge-SET-current' in- 
stead of the usual direct SQUID scheme 'magnetic flux- 
current'. We have analyzed the charge sensitivity as a 
function of magnetic flux, gate voltage, Coulomb inter- 
action, dot asymmetry, and temperature. The sensitiv- 
ity of our setup can be further increased by adding an 
electromechanical element [16]: Applying a large elec- 
tric field to the charged nanowire, the change in charge 
will lead to a mechanical shift of the wire. This shift can 
then be detected due to the change in the capacitance 
of the compound setup as in Ref [16]. In the present 
work, we have concentrated on a single-channel wire in 
order to demonstrate the effect; the case of an n-channel 
wire (n = 2 or n > 2) can be analyzed using the same 
technique and we plan to study this case in the near 
future. 
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